Master singular behavior for the Sugden factor of the one-component fluids near their 
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We present the master (i.e. unique) behavior of the squared capillary length - so called the 
Sudgen factor-, as a function of the temperature-like field along the critical isochore, asymptotically 
close to the gas-liquid critical point of twenty (one component) fluids. This master behavior is 
obtained using the scale dilatation of the relevant physical fields of the one-component fiuids. The 
scale dilatation introduces the fluid-dependent scale factors in a manner analog with the linear 
relations between physical fields and scaling fields needed by the renormalization theory applied to 
the Ising-like universality class. The master behavior for the Sudgen factor satisfies hyperscaling 
and can be asymptotically fitted by the leading terms of the theoretical crossover functions for 
the correlation length and the susceptibility in the homogeneous domain recently obtained from 
massive renormalization in field theory. In the absence of corresponding estimation of the theoretical 
crossover functions for the interfacial tension, we define the range of the temperature-like field where 
the master leading power law can be practically used to predict the singular behavior of the Sudgen 
factor in conformity with the theoretical description provided by the massive renormalization scheme 
within the extended asymptotic domain of the one-component fiuid "subclass". 

PACS numbers: 64.60.-i, 05.70.Jk, 64.70.Fx 



1. INTRODUCTION 

The knowledge of interfacial properties [l| for a non- 
homogeneous fluid of coexisting vapor and liquid at equi- 
librium is of prime importance for many engineering ap- 
plications and process simulations. Moreover, accurate 
predictions of these interfacial properties are essential to 
gain confidence in modeling underground geological fluid 
flows in porous media, oil recovery, gas storage in geo- 
logical formations, pool boiling phenomena, microfluidic 
devices based on wetting phenomena, etc. 

A large number of different forms of related phe- 
nomenological laws, the so-called ancillary equations, are 
reported in the literature to calculate interfacial prop- 
erties along the vapor-liquid equilibrium (VLE) line [2]. 
These relations complement the complex multiparameter 
equations of state (EOS's) which have been developed to 
accurately fit the thermodynamic properties measured 
in the homogeneous domain. Such a phenomenologi- 
cal approach to estimate fluid properties is commonly 
based on the multiparameter corresponding-states prin- 
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ciple 0, i, i|. In the following we call n-CSEOS such 
an EOS which contains n > 2 system-dependent param- 
eters. The main reason for the power of such a phe- 
nomenological approach is related to the fact that the 
two-parameter corresponding-states (2-CS) principle can 
be applied to any polynomial EOS which has a liquid- 
vapor critical point However, in spite of increasing 
the number of fluid-dependent parameters, the common 
calculation of interfacial properties from ancillary equa- 
tions and n-CSEOS, is not only mathematically complex, 
but is also unable to account for: 

1) the molecular fluid complexity Q, especially the 
non-spherical symmetry of molecules and the quantum 
behavior of light fluids Q ; 

2) the asymptotic scaling of the critical phenomena 
close to the gas- liquid criticalpoint [6j, especially the 
non-analytic Ising-like nature 0, [3| of the critical expo- 
nent 

Among these interfacial properties, the capillary length 
ica, or more precisely the squared capillary length {tca)'^ 
also called the Sugden factor and noted Sg in the 
following, plays a special role on Earth's gravity envi- 
ronment (recalled here by the subscript g). The Sugden 
factor reflects the balance between interfacial and volu- 
mic forces which defines the shape and the position of the 
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interface in equilibrium when subject to the gravity field 
of constant acceleration g. In the case of perfect Hquid 
wetting, S'g is then related to the surface tension a and 
the density difference Ap^v = PL — Pv between coexist- 
ing liquid (density pl) and vapor (density pv) phases by 
the equation 



2a 



- gAp 



LV 



(1) 



where a is the surface tension. Therefore, the knowledge 
of the Sugden factor is an important challenge to provide 
better control on non homogeneous fluid properties. 

In addition, as clearly documented two decades ago 
[ill . [l3, [l3| , the temperature dependence of Sg , along a 
large temperature range of the VLE line of all investi- 
gated one-component fluids [III, [H, [li, [H, [H, [13, [H, 
[lol . Ilol . [21] . [2^ , shows a pure power law behavior which is 
applicable over an appreciably larger temperature range 
[see below Eq. |(2]) and the related dicussion of the Fig. 
[1^]. Such a weak temperature dependence of the effective 
exponent at small but finite distance of the critical point 
was partly well-understood to be related to the smallest 
value (~ 0.51 [§|, see below) of the confluent exponents 
which govern the corrections to asymptotic scaling of crit- 
ical phenemonena [2l|- However the theoretical reason to 
observe a near zero- value of the amplitude contribution 
of the confluent corrections for the Sugden factor case 
remains unclear, specially in the absence of estimation of 
the crossover behavior of the surface tension. 

Indeed, the significant theoretical improvements to ac- 
count for classical-to- critical crossover Q, specially in 
the one-component fluids 0, provide the most power- 
ful tools available today to analyze accurately interfa- 
cial properties in large temperature ranges of the non- 
homogeneous domain. For example in the present work, 
using the crossover functions recently derived 
from the massive renormalization scheme [IE, [13 
our main objective is to accurately estimate this lead- 
in g asy mptotic behavior of Sg from scalin g ar guements 




0^,|31|,|3^ and available MR description [3^,|34,|35j of 
the master singular behavior of the one-component fluid 
"subclass". Such a description is based on the formal 
analogy between the scale dilatation of the physical fleld 
variables proposed by Garrabos [3B,[33,[3l| and the Hnear 
relation between the physical fields and the scaHng fields 
needed by the renormalization theory [s^. The major 
advantage of this scale dilatation method is to estimate 
the universal behavior of any one-component fiuids with- 
out adjustable parameters, by using only the four critical 
coordinates of its liquid-vapor critical point (excluding 
here quantum fluids jsH] to simplify the presentation of 
the scale dilatation method). 

The paper is organized as follows. Section 2 demon- 
strates the master singular behavior of Sg observed from 
the scale dilatation method. The corresponding Ising- 
like asymptotic description of Sg based on hyperscaling 
0,[33,[3l|,l32| and MR description [ii,^,^! of the criti- 
cal crossover is reported in Section 3. The master leading 



terms of the MR crossover functions for the correlation 
length and the susceptibility in the homogeneous domain 
[24 . , are used to demonstrate that the fit of the mas- 
ter behavior observed in the (nonhomogeneous) extended 
asymptotic domain can be made with a theoretical preci- 
sion of the same order of magnitude than the experimen- 
tal one. The discussion given in Section 4 shows the main 
points to be considered for a classical-to-critical crossover 
description of the interfacial properties at finite temper- 
ature distance to the critical temperature. Specifically, 
we estimate precisely the temperature-Hke range where 
this theoretical treatment becomes unappropriate to rep- 
resent the increasing non critical microscopic difference 
between gas and Hquid approaching the triple point tem- 
perature. Conclusion is given in Section 5. 



MASTER SINGULAR BEHAVIOR OF THE 
SUGDEN FACTOR 



2.1. The data sources 

The Sugden factor measurements Sg (|Ar|) as a func- 
tion of the temperature distance AT = T — Tc to the 
critical point in the nonhomogeneous range T <Tc, have 
been published and analyzed for several one-component 

fluids [13, [H, lii, [H, M, [13, E, E, [13, [HI, [13,1131 . t 

(Tc) is the temperature (critical temperature). Sg (|AT|) 
data are generally obtained along the critical isochore 
/9 = Pc in a flnite temperature range bounded by max 

T, where Tc,exp is the 



and min values of I AT I — T, 



c,exp 



measured (or estimated) critical temperature in the ex- 
periments [p (pc) is the density (critical density)]. The 
relative precision clamed by the authors is generally lower 
than 10%. For most fluids, Sn was fltted using the effec- 
tive power law of equation [J [HI, [l3| 



Sg — So^e 



(2) 



where the dimensionless temperature distance |At*| to 
the critical point was defined by 



|At*| = 



lATl 



T 



T, 



c,exp 



T, 



(3) 



c^exp 



In Eq. |[2]), the free amplitude So,e was a fiuid-dependent 
quantity related to the effective value tpe — 0.91 — 0.97 
of the free (or fixed) exponent ipe considered as an ad- 
justable parameter when measurements were performed 
in a restricted temperature range at finite distance to 
Tc,exp- The corresponding results ipe;So^e [with free (or 
fixed) value of ipe] for each selected fluid are summarized 
in columns 3 and 4 of Table [J (references are given in 
column 2). However admitting now that At* is the rel- 
evant physical fleld |[30j to describe the singular scaling 
behavior of the thermodynamic fluid properties in the ho- 
mogeneous or nonhomogenous domains along the critical 
isochore, the three main "critical phenomena" features of 
these fltting analyses are: 
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i) the correlation between the effective values of 5*0,6 
and ipe is highly dependent on T^.exp and on the (min 
and max) values of the temperature range covered by 
the fit "close" to the critical point; especially when the 
local values of are estimated in the temperature range 
lower than |Ar*| < 0.05, the common averaged value 
(f>e = 0.935 ± 0.015, equal to the asymptotic universal 
value tp = 2v — (3 = 0.935 ± 0.015 obtained from the 
present theoretical estimation of the critical exponents 
V = and /? = |i| [see below Eq. {H)], appears consistent 
with the Sg data, whatever the one-component fluid. 

ii) accordingly, the temperature dependence of the ef- 
fective exponent over a larger temperature range is very 
small (then with a sign not unambiguously deflned for the 
small amplitude of the leading confluent term) , whatever 
the one-component fluid or the extension of the temper- 
ature range of the flt; 

iii) the measured value of the effective exponent is 
never equal to the mean-fleld value V'mf — 1, whatever 
the one-component fluid and even the large temperature 
range of the flt; 

As a matter of fact, it is well-established now [S] that 
the range of validity of the asymptotic scaling form of Eq. 
([2]) is, strictly restricted to the asymptotic approach of 
the liquid-gas critical point (CP), when a oc \At*\'^ and 
ApLv oc |Ar*|'' simultaneously go to zero for At* — > 
in Eq. JT]). = (d- l)r^ « 1.261 and P « 0.326 are 
the universal values j§] of the critical exponents related 
to a and Aolv, respectively, and v ~ 0.63 is the univer- 
sal value [9(1 of the critical exponent for the correlation 
length, with ^ oc |Ar*|~'^ (see also Q for details of nota- 
tions and deflnitions). However, at small but flnite |At*|, 
any pure power law like Eq. ^ must be modifled to ac- 
count for confluent corrections to scaling which can be 
represented by the Wegner-Hke expansion with the 
universal features of uniaxial 3D Ising like systems 
Then, the asymptotic singular decrease of Sg must be 
fltted by the following equation 



Sg = So I At* 



l + J2s,\Ar*\ 



i=l 



(4) 



where A « 0.51 is the universal value [§| of the critical 
exponent which characterizes the leading family of the 
confluent corrections to scaling. The amplitudes 5*0, ^i, 
... Si, etc., are fluid-dependent quantities. Equation (0]) 
means that the critical exponent 



ip^(l)-l3^{d^l)iy-f3 



(5) 



only takes its universal value ip ~ 0.935 asymptotically 
when At* — > 0. Therefore, the weak temperature depen- 
dence of the effective exponent at flnite value of At* , flrst 
shows low rate of convergence of the Wegner expansion. 
Moreover, in the fltting of the experimental Sg data, if 
the contribution of the confluent correction terms is neg- 



ligible, then y^5i 



At 



in Eq. g]). 



To illustrate this behavior, the Sugden factor Sg = 
i^Ca)^ (expressed in w?) can be divided by (Tc — TY 
(expressed in if"^, with p — 0.935) [l3|. In Figure[lK (log- 



log scale) , this convenient scaled form 



is shown as a function of the temperature distance Tc — T 
\K] for eighteen one-component fluids. Each curve has a 
relative temperature extension corresponding to the ex- 
perimental temperature range (including for some fluids 
measurements until their triple point temperature Ttp). 
The use of such dimensional quantities makes the order 
of magnitude of the leading amplitude contribution (i.e. 
Sg {T — Tc + IK) ~ So (Tc) '^) of each fluid clearly dis- 
tinguishable, while the quasi-horizontal line whatever the 
fluid (except the water case which needs a special atten- 
tion given in § 4.3) shows that the confluent contribution 



is negligible (i.e. y^S^ (Tc) 



1 



(T, - Tf'^ 



0). Sg val- 



T — 1 K cover one decade: from 1.5 10 ^ 

3.84 mm^], to 



ues at Tc 

for sulfurhexafluoride [with So (SFq) 
1.5 lO-'^m^ for methane [5*0 {CH4) = 13.80 mm^]. 

We have noted that some of the data reported in Fig. 
[IK have been measured in a large temperature range of 
the coexisting VLE line, including measurements close 
to Ttp. To separate the asymptotic critical range from 
the triple point location along the temperature axis, we 
have marked by vertical arrows the temperature distance 
where T = 0.7 Tc (i.e. the temperature distance where 
the practical fluid-dependent acentric factor [4l| is de- 
flned in the nonhomogeneous domain) . The "non-critical" 
temperature range between 0.3 Tc < T^ — T < Tc — Ttp 
(the right hand side of the corresponding arrows in Fig. 
dK) is considered to be far away from the critical point. 
Close to the critical point, the practical temperature 
range where the Wegner-like expansion flt the singu- 
lar behavior does not usually exceed a few percent in 
I At* I [42[. In a similar arbitrary manner, we have repre- 
sented by vertical arrows the temperature distance where 
T — 0.99 Tc, to make clear the "critical" temperature 
range Tc — T < 0.01 Tc (the left hand side of the cor- 
responding arrows in Fig. [l^), where the use of Eq. (0]) 
has a theoretical justiflcation, as discussed below in § 
3. To introduce the main physical parameters needed 
for accurate description of the singular behavior of Sg in 
this asymptotic temperature range, the next subsection 
presents the application of the scale dilatation method 
p^ . [37! leading to deflne the dimensionless form (noted 
Sg) and the renormalized form (noted S*, ) of the Sugden 
factor, with two objectives: 

1) to show that any modehng based on the 2-CS prin- 
ciple is inaccurate to describe the fluid dependence of 
(dimensionless) S* [see below Eq. (fT2|) ] as a function of 
the (dimensionless) temperature fleld At*; 

2) to unambiguously show the master singular behavior 
of (renormalized) S*, [see below Eq. ifTS)) ] as a function of 
the (renormalized) temperature-like fleld, noted T* [see 
below Eq. ^] . 
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Fluid 


Ref 




So,e 


So, if Ref 


Zs,ip 


&Zs,ip 
% 


Ar 


[17J 
[12] 


0.940 
0.913 


4.13 
3.78 


4.036 [13| 
4.18 This work 


2.423 
2.510 


-1.9 
+1.6 


Xe 


[18] 


0.942 


3.05 


2.953 [1^1 


2.668 


+8.0 


N2 


[17J 
[12] 


0.930 
0.926 


5.46 
5.10 


5.59 This work 
5.32 This work 


2.446 
2.328 


-1.0 
-5.7 


O2 


[12] 


0.909 


4.85 


5.47 This work 


2.563 


+4.6 


CO2 


[19] 


0.933 
0.920 


9.47 
8.40 


9.52 \13} 
9.0 This work 


2.55 
2.411 


+3.3 
-2.4 


SFe 




0.943 


3.931 


3.84 [13] 


2.46 


-0.4 


CCI3F 




0.928 


6.234 


6.44 This work 


2.470 


0.0 


CCI2F2 




0.936 


5.615 


5.589 This work 


2.476 


+0.3 


CCIF3 


[19] 


0.972 
0.9379 


5.33 
4.847 


4.5 This work 
4.783 This work 


2.268 
2.411 


-8.1 
-2.4 


CBrF'i 




0.938 


3.879 


3.826 This work 


2.374 


-3.1 


CHCIF2 


M 


0.937 


6.859 


6.796 This work 


2.323 


-6.0 


C2H4 


M 






13.90 [13] 


2.480 


+0.4 




m 


0.933 


13.6 


13.73 This work 


2.382 


-3.6 


C2HQ 








14.42 [13] 


2.437 


-1.3 


i — CiH\o 


[15] 






12.71 [13] 


2.392 


-3.2 


n - C5H12 


[40] 


0.935 


12.916 




2.488 


+0.7 


n - CaH-iA 


[40] 


0.935 


12.753 




2.445 


-1.0 


n - C7H16 


[40] 


0.935 


12.520 




2.552 


+3.4 


n - CsHis 


[40] 


0.935 


12.217 




2.520 


+2.0 


H2O 


[22] 
[21] 


0.91 


33.2 


34.72 [13| 
37.25 This work 


2.262 
2.427 


-8.4 
-1.7 












2.4530 
2.47 


±3.1 



Table I: Effective leading amplitude So,e [see Eq. ^] of the Sugden factor Sg = {Ica)'^ {£ca is the capillary length) (column 
2) for several one-component fluids (colum 1); Calculated values of the master amplitude-Es.o (column 3) of the renormalized 

Sugden factor S** [see Eq. lfl8|l ]. using Eq. (Ell); The residual SZs = 100 x (^^-^f^ - l) (expressed in %) from the value 

Zs = 2.47 estimated from Eq. I|43p is given in column 4; (for data sources and the selected fitting results see the references 
given in the last column). 



2.2. The scale dilatation method to observe the 
master singular behavior 

The following analysis of the Sugden factor from the 
scale dilatation method is similar to the one of the corre- 
lation length given in Ref. [33|. We recall only the main 
features (ignoring the quantum contributions at T = Tc 
[3^). The input data are the four critical coordinates 



dp 
df 



(6) 



which localize the liquid gas critical point on the phase 
surface of equation <I>2_ {p,Vp,T) = for each fluid par- 
ticle of mass mp p is the pressure, Vp is the particle 
volume, and ap {T, Vp) is the Helmholtz energy per parti- 
cle. The subscript p refers to a particle quantity and all 



the definitions and notations related to Eq. ^ are given 
in [sB, [13, [Hi . The critical data related to the fiuids se- 
lected in Table [J are reported in Table [ill We note that 
Tc values of Table (llj which were obtained from the ther- 
modynamic analysis of the phase surface, can be slightly 
different from Tc,exp values given in the experiments ref- 
ered in Table HI Also pc — values from Table Hi] can 
be slightly different from the experimental critical den- 
sity values reported in these experiments. 

In combining Q™", the Boltzmann constant /cb, and 
d = 3, Eq. |[6]) can be written in a more convenient form, 
such that 

= {(/3c)"\ac,^c,yc} (7) 
which introduces the following four scale factors given by, 

(/3c)"' = fcsT, (8) 
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Fluid 


mp 


Tc 
(K) 


Vp,c 


Pc 

(MPa) 


7c 

{MPaK-^) 


(lo-^V) 


(nm) 


z,. 




Ar 


6.634 


150.725 


0.12388 


4.865 


0.191 


2.08099 


0.753463 


0.2896 


4.32882 


Xe 


21.803 


289.733 


0.19589 


5.84 


0.1182 


4.0003 


0.881508 


0.28601 


4.85434 


N2 


4.652 


126.214 


0.14814 


3.398 


0.1715 


1.74258 


0.80043 


0.288875 


5.37014 


O2 


5.314 


154.580 


0.12187 


5.043 


0.1953 


2.13421 


0.750786 


0.287972 


4.98641 


CO2 


7.308 


304.137 


0.15622 


7.3753 


0.170 


4.19907 


0.82882 


0.27438 


6.0104 


SFe 


24.255 


318.735 


0.32769 


3.754 


0.0835 


4.40062 


1.0544 


0.27954 


6.0896 


CCI3F 


22.810 


471.110 


0.41174 


4.4076 


0.0655 


6.50438 


1.13850 


0.27901 


6.00530 


CCI2F2 


20.078 


384.930 


0.35562 


4.1249 


0.0745 


5.31454 


1.08814 


0.27602 


5.95186 


CClF-i 


17.348 


301.88 


0.29807 


3.877 


0.0910 


4.16791 


1.02441 


0.27727 


6.08861 


CBrFz 


24.727 


340.19 


0.33191 


3.956 


0.0810 


4.69683 


1.05889 


0.27956 


5.96985 


CHCIF2 


14.359 


369.30 


0.27454 


4.990 


0.0965 


5.09874 


1.00721 


0.26869 


6.14259 


C2H4 


4.658 


282.345 


0.21667 


5.042 


0.11337 


3.89820 


0.91781 


0.28131 


5.34856 


CH4 


2.664 


190.564 


0.16361 


4.5992 


0.14442 


2.63102 


0.830133 


0.28679 


4.9838 


C2He 


4.993 


305.322 


0.24171 


4.872 


0.10304 


4.21554 


0.95290 


0.27935 


5.45505 


i — C4H10 


9.652 


407.844 


0.43020 


3.629 


0.0610 


5.63084 


1.15770 


0.27726 


5.93407 


n - C5H-12 


11.981 


469.70 


0.521785 


3.3665 


0.0511 


6.48491 


1.24425 


0.270875 


6.12956 


n - CgHu 


14.310 


507.49 


0.61138 


3.0181 


0.043658* 


7.00666 


1.3186 


0.26667 


6.30719 


n - CrHie 


16.6386 


540.13 


0.7168 


2.727 


0.038068* 


7.45731 


1.3983 


0.26218 


6.64356 


n - CsHis 


18.9683 


568.88 


0.81839 


2.486 


0.033768* 


7.85424 


1.46746 


0.258978 


6.82776 


H2O 


2.991 


647.067 


0.09268 


22.046 


0.275 


8.93373 


0.740 


0.229 


7.071 



Table II: Set of critical parameters [see Eqs. (HJ and ^\ for twenty one-component fluids of particle mass mp [43] selected in 
the present work (see Fig. 1). 



Zr 



ksTc 



Pc 



(9) 



(10) 



(11) 



Equation |[7]) involves one energy scale unit (/3c) , one 
length scale unit ac, and two dimensionless scale factors 
Zc and Yc characterizing two preferred directions to cross 
the critical point along the critical isotherm and the crit- 
ical isochore, respectively, etc, which does not depend of 
the size L ^ iV) '' of the container, has a clear physical 
meaning as length unit js^: it represents the spatial ex- 
tent of the short-ranged (Lennard- Jones like) molecular 
interaction [31, which allows us to define Vc.i = -^^^ 
as the volume of the microscopic critical interaction cell 
(CIC) of each fiuid. Zc is the usual critical compression 
factor. Furthermore, (Zc) = ricVcj is the number of 
particles that fills Vcj. Then the minimal set of data in 
Eq. ((7| is related to the thermodynamic properties of 

the critical interaction cell of size Uc = {vcj)^ (45l |. 

We recall that the critical compression factor Z^ and 
the critical Riedel factor aB.,c [related to Yc by aR,c = 



Fc + 1] , are among the basic parameters used to develop 
4-CSEOS's for engineering fiuid modeling |47i] . 

The characteristic units {(}c)~^ and ac are the param- 
eters needed to provide a dimensionless analysis of the 
fiuid properties, leading to their "classical" description 
based on the two-parameter corresponding state (2-CS) 
description. Obviously, the dimensionless form of the 
Sugden factor is given by 



Sn 



(12) 



Figure [Tb (log-log scale; color online) represents the con- 

fiuent behavior of the rescaled dimensionless quantity 

s* 

| y, as a function of the dimensionless temperature 
distance At*. Figure [TJd complements Fig. 3 initially 
published by Moldover in Ref. after normalization 

of the vertical axis by (ac)^- Figure [TJd illustrates the re- 
sults of any classical two-parameter corresponding state 
theory (here the two characteristic parameters are (/3c) 
and Qfc). Figure [TJd, shows the failure of the 2-CS prin- 
ciple in terms of molecular fiuid complexity since, from 
xenon to water, the dimensionless Sugden factor covers 
one order of magnitude at the same reduced tempera- 
ture distance to the critical point. Moreover, in terms 
of classical critical phenomena, using Eq. |(T]) where 
ApLv oc |At*|''*^^ and a cx |Ar*|"*"-^^ with /3mf = 5 
and (pMF = I (3) obtain the mean field exponent 
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'^MF = 1- This mean-field value associated to the clas- 
sical behavior of the correlation length (with exponent 
vmf = \) expected from Van der Waals-Hke theories 
[i^ . [49| . is unable to describe the experimental results, 
even at large temperature distance, as shown by the sig- 
nificantly positive slope fMP — f — 0.065 reported in Fig. 
[Hd. In addition, the scaling law {d — \)v = 4)^ that ex- 
plicitly involve d, is not correct for mean-field exponents 
in three dimension. We will turn back on the mean- 
field theories in § 4.2 when we will discuss the related 
critical-to-critical crossover description of the interfacial 
properties. 

In the next step, the dimensionless scale factors Yc and 
Zc are introduced throughout the scale dilatation method 
[131 ■ Typically, the scale dilatation of the dimensionless 
temperature distance. 



Ar* = ksPc (T - Tc) 
leads to the renormalized thermal field, 

T* = KAt* 



(13) 



(14) 



The scale dilatation of the dimensionless order parameter 
density 

Am* = (ac)''(n-nc) = (Zc)"'Ap (15) 
leads to the renormalized order parameter density 

M* = {ZcY Am* = (Zc)^ Ap (16) 

In addition, the renormalized form £*=£* = — of the 
correlation length ^ [3H|, leads to the renormalized form, 
noted E*, of the surface tension a such that [50|| 



E* = cr* = af3c (ftc) 



d-l 



(17) 



Taking into account Eqs. (U) and (fT2l) . the renormalized 
Sugden factor S** reads as follows |i50|] 



5* =g*(Z,) H^ca) 



d-l 



(18) 



with g* = rupPcCtcg- Therefore, after application of the 
scale dilatation method, the renormalized form of Eq. (U) 
is given by 



2.3. Predictive power of the scale dilatation 
method within the Ising-like preasymptotic domain 

As initially shown in Ref. we can expect to fit the 
master singular behavior of 5*. observed asymptotically 
close to the critical temperature by a restricted (two- 
term) Wegner-like expansion given by 



Zs IT* 



1 



Zl\T* 



(21) 



where 4> ~ 0.935 and A « 0.51 are the universal crit- 
ical exponents while Zg and Zg are the master (i.e. 
unique) leading and confiuent amplitudes, respectively, 
for all one-components fiuids. By term to term compari- 
son of Eqs. (0]) and l(2T1) using Eqs. l(20l) . we obtain the 
following relations 



2s = g' 



(ttc)^ (Zc 



So 



Zl^{Y,)-^S, 



(22) 



(23) 



which shows the unequivocal link between master ampli- 
tudes and system-dependent amplitudes, through Q™'" 
[(Eq. m\- 

In other words, only when the fiuid-dependent set 
Q™'" and the master amplitudes Zg and Zg are 
known, the restricted Wegner-Hke expansion [Eq. (|4]) 
with i < 1] o{ Sg can be determined for any one- 
component fluid by inverting Eqs. l|22p and i|23p . such that 

^0 = (acf-^ {g*)-^ {Z,Y {Y^yZs and 5i = {Y^f 2%. 
Then, the master values of Zs and Z\ conform to the 
universal features calculated for the Ising-like universal- 
ity class (i.e., some combinations and ratios of Zs and 
Zg take universal values, in agreement with the two- 
scale-factor universality). We will detail this point in § 
4. Before, in the next Section, the scale transformations 
of Eq. (|20| are reported in conformity with the asymp- 
totic linearization [3^ of the two relevant fields needed 
by the renormalization group theory. That leads indeed 
to the correct account for universal features estimated in 
the preasymptotic domain and the accurate determina- 
tion of Zs using the present theoretical status provided 
by the MR scheme 



M 



(19) 



LV 



As expected [3d[, the collapse on the master curve ob- 
tained from the scale transformations 

Ar* ^ T* = YcAt* 

S; ^ S;. = g* {Z,)-i S; (T1 case) (20) 
-^5g*. X IT*^^ {T2case) 

is shown in Fig. [1]: {Tlcase) and [TJi {T2case), inde- 
pendently of any theoretical form used to represent this 
master behavior. Now the scatter of the collapsed data 
corresponds to the estimated precision (10%) for the Sug- 
den factor of each fluid. 



3. ISING-LIKE CROSSOVER FUNCTIONS FOR 
THE SUGDEN FACTOR 

To our knowledge, the theoretical function giving the 
classical-to-critical crossover of the interfacial tension 
a ( Ar* ) is not available from the MR scheme, while the 
one of the coexisting density ApLv (At*) [23| remains 
affected by a large uncertainty on the value of the first 
confluent amplitude. Therefore, using either Eq. ^ for 
physical properties or Eq. I|19|) for renormalized proper- 
ties, the related crossover functions of the physical and 
renormalized Sugden factor remain undetermined. Es- 
pecially the value of Zg [Si, respectively] in Eq. l(2T|) 
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Figure 1: (Color online) a) Singular behavior (log-log scale) of ^.^ ^y^y (expressed in K with </5 = 0.935), as a function 
of the temperature distance Tc — T, for nonhomogeneous one-component fluids (see Tables I and II). For each fluid, the arrows 
indicate the arbitrary temperature distances of Tc — T = O.OlTc (left) and Tc — T = O.STc (rigth); Inserted Table gives the color 
indexation for each fluid (see Table I and text for details); b) Log-Log scale of the respective dimensionless variables using as 
a length unit and (/3c)~^ as a energy unit. Each distinguishable curve illustrates the failure of the classical corresponding state 
scheme. The expected slopes of a classical power law with mean field exponent ipMF are given by the directions labelled MF 
(see text); c) Master singular behavior (log-log scale) of the renormalized Sugden factor 5*. — g* (Zc)'^ S* (with S* = j-§yi 
and g* = mpfScCtcg), as a function of the renormalized thermal field T* = Yc |Ar*| [see text and Eq. I|20p ]: d) Master "confiuent" 
behavior of the rescaled quantity yf^ryp, as a function of the renormalized thermal field T* = Yc | Ar*|, [see text and Eq. lf20|l [. 
In c) and d): The curve labelled MR^ corresponds to the Eq. l[40|l using crossover Eqs. l[32|l and ifSSjl for the correlation 
length and the susceptibility, respectively, in the homogeneous domain; d'p^]^ [Eq. dM))] and [see Eq. jHH correspond 

to the extension of the preasymptotic domain and the extended asymptotic domain, respectively; The graduation of the upper 

horizontal axis gives — — calculated from master crossover Eq. l[32|l : At large values of the renormalized thermal field 

(T* >o) 

(T* > 0.2) which correspond to — — J; li the restricted range labelled Vec indicates the effective crossover for the 

exponent (3, while the vertical index labelled m gives a rough estimate of the microscopic limit — ^-^ ^ ~ \ where the 
correlation length ^ in the nonhomogeneous domain will reach the order of magnitude of the two particle equilibrium position 
re (with re ^ cr, where a is the size of the particle); See inserted Table and legend of a) for fiuid color indexation; See legend 
of b) for arrows in lower horizontal axis. 



8 



[Eq. ([4]), respectively] cannot be estimated from theo- 
retical prediction of the universal value of the confluent 
amplitude ratios related to the lowest confluent exponent 
A (see also our discussion in Section 4.1). However, hy- 
perscaling related to the two-scale-factor universality of 
the asymptotic Ising-like description provides unambigu- 
ous determination of the values of [Sq, respectively] 
in Eq. l(2T|) [Eq. respectively]. This determination 
is presented below using the master forms of Ising-like 
crossover functions obtained from the massive renormal- 
ization (MR) scheme. 

However, we note that a form equivalent to Eq. i(2T|) 
was also recovered in the crossover approach of Belyakov 
et al [s^l, who uses adjustable parameters as scale factors 
of the physical variables. The solution was obtained on 
the basis of the e-expansion in first order e and was not 
considered here due to the arbitrary of the phenomeno- 
logical adjustement to provide the crossover to a classical 
behavior. 



3.1. Asymptotic hyperscaling description of the 
Sugden factor 

It is well-established experimentally [l^ . [T3 | and the- 
oretically P, m, [H, HBl, that the asymptotic Hmit for 
At* — > of the product of the interfacial tension by 
the squared correlation length takes a universal value, 
noted i?^^, for the Ising-like universality class. This re- 
sult proceeds from the Widom's scaling law between the 
corresponding critical exponents v and /i given by 



(24) 



with d = 3 in our present study. Therefore, we can in- 
troduce R^c as follows 



i?± =/3,xlim a(|AT*|) X K(Ar*)] 



Ar*^0± 



(25) 



singular properties 



iPc (p. 

lim 



KT(Ar-)x[C(Ar*)]- 
[Ap^v(IAr-l)]- 



(27) 



*0± 



Equation l|27p relates the singular behaviors of the cor- 
relation length ^ (Ar*) [with critical exponent u and 
leading amplitude ^q] in the homogeneous domain, the 
isothermal compressibility kt (At*) [with critical expo- 



nent 7 and leading amplitude Tq 



-] in the ho- 



mogeneous domain and the order parameter density 
ApLv [with critical exponent f3 and leading am- 

plitude Bq = 2pcB] in the non-homogeneous domain. 

Using Eqs. |(T|), l(25l) . and l(27|) to eliminate both prop- 
erties (T(|Ar*|) and ApLv {\At*\) of nonhomogeneous 
domain, we obtain the following asymptotic equation 



lim [S, 



lim 



2(/3c)'Pcff 



[KT(Ar*)xC(Ar*)]2 



Ar*^0+ 



(28) 

which relates the asymtotical singular behavior of the 
Sugden factor in the nonhomogeneous domain to the ones 
of Kt (At*) and ^ (At*) in the homogeneous domain. 
The corresponding scaling law reads 



d .\ 7 



(29) 



The scaling laws given by Eqs. 1^^, ([26]), and where 
explicit reference to the space dimension is needed to 
connect correlation exponents and thermodynamic expo- 
nents, are characteristic of hyperscahng and reflect the 
universal features related to the two-scale-factor univer- 
sality, which do not depend on the (homogeneous or non- 
homogeneous) domain (see also Ref. j56.]). 



where the superscript ± refers to the singular behavior 
of ^ above (-I-) or below (— ) Tc. As a matter of fact, ac- 
counting for the universal ratio ||Ar'-<o) =1-96 for the 
Ising-like universality class jsj, the amplitude combina- 
tion i?Jr^ — (1.96)'^ ^ shows that an interfacial prop- 
erty (here a cx [At*!"^) in the non-homogeneous domain 
(At* < 0) is related in a universal manner to the corre- 
lation length in the homogeneous domain (At* > 0). 
Considering the scaHng law 

dv = -f + 2P (26) 

it is also well-estabHshed that the amplitude combina- 

tion ( ^ j -gT , noted , j , (using customary nota- 

tions [7|) corresponds to the universal value of the asymp- 
totic limit for At* — > of the following combination of 



3.2. The master crossover of the one-component 
fluid subclass 

We are now able to construct one pseudo-crossover 
function based on Eq. I|28p . This pseudo-crossover 
function for the Sugden factor accounts exactly for the 
asymptotic two-scale factor universality but agrees only 
qualitatively with the one-parameter Ising-like critical 
crossover description at finite distance to CP. As a matter 
of fact, accurate expressions of the complete classical- 
to-critical crossover were recently proposed by Bagnuls 
and Bervillier ^2j| and written in appropriate Ising-Hke 
asymptotic forms by Garrabos and Bervillier [2^ to ac- 
count for error-bars associated with the estimations of the 
universal exponents near the non-Gaussian fixed point. 
Moreover, introducing only three characteristic numbers, 
L{i/}, e^^f\ and (see Ref. Q for details). 
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these crossover functions can be easily modified to accu- 
rately describe the master singular behavior of the one- 
component fluid subclass. In this master description, two 
leading amplitudes Z+ , , and one confluent amplitude 

among Z^''^ and Z^''^ , can be selected as characteristic 
parameters of the Ising-like universal features observed 
in the Ising-Hke preasymptotic domain. Z+, Z^ , ^x^^ 

and Z^''^ are associated to the asymptotic crossover be- 
havior of the correlation length and the susceptibility 
in the homogeneous domain. We recall that the corre- 
sponding master crossover functions are asymptotically 
approximated by the restricted (two terms) Wegner like 
expansions given by the respective equations 



t (T* 



K* (T* 



z+(r*)- 



Z+ [T* 



1 



(T* 



(30) 



(31) 



where Z+ = 0.119, Z+ = 0.570, 2^^+ = 0.555 and 
Z'^''^ = 0.377 are the constant values of the master (i.e. 
fluid independent) amplitudes, with the universal ratio 

-A-Tf = 0.68 [9]. Accordingly, the modifled crossover 



functions are given by the following equations 



i=3 



(32) 



j=3 



(33) 



with 



and 



D{t) 



l + S2Vi 



t = Ql^f^T* 



(34) 



(35) 



All the critical exponents {ly, 7, A, Amf) and the con- 
stants (Z^, Z+, X^^i, Y^^i, X-^^i, Y^^i, S2) of the initial 
crossover functions defined in Ref. [IBl are reported in 
Tableini Furthermore, in Eqs. ^ and the prefac- 

tors and ZJ. relate the asymptotic master behav- 

ior given by Eqs. ([30l) and ^Bl, respectively, and satisfy 
to unequivocal estimations from the three characteristic 
numbers L,^^^^, 8'^"'^^^, and vj/'fi/^ of the one-component 
fiuid subclass fsal, such that, 



Z 



{1/} _ 



Z+Z+ 



ell/} 



(36) 



Z 



{1/} 



Z+Z+ (e{i/>) 



(37) 



The scale factor Q^^f^ is defined from the following ratios 
of the confiuent amplitudes 



Q{lf} 




(38) 



where Z, 



-ES^x,»^x.«> with = 0.68 125)]. All the 

values of these master constants are shown in Table IIVI 

1 1 f I 

We also note that the master prefactors ZJ and 

{1/} 



and 
0.68 125)]- 



z 



, as all the other prefactors which modify the initial 



crossover functions to account for master behavior of the 
renormalized properties of the one-compenent fiuid sub- 
class, take the same value above and below the critical 
temperature, while only two of them are characteristic 
of this subclass. In addition, the single master crossover 
parameter Q^^^^ is the same for any property along the 
critical isochore, above and below the critical tempera- 
ture. As demonstrated in Refs. [H, [S^, it is possible to 
define unambiguously the extension T* < CpJ^^ of the 
preasymptotic domain where each master crossover func- 
tion can be approximated by its restricted (two-term) 
expansion. Using e{i^> (see Table Hv]) we obtain 



n-* < A^f^ — '-PAD 
^ ^ '-PAD 



10 



-3 



where C 



I sing 
PAD 



{S2 



is defined in Ref. |2£ 



10" 



(39) 



After appropriate rescaling of the master form of each 
property included in Eq. (|28l) . we define the following 
master quantity 



SiT* 



R 



>t* (T*) t (T*) 



(40) 



where the correlation length and the susceptibility are 
given by Eqs. {M]) and ([Ml), respectively. SiT*) 
[Eq. l(40|) ] is the pseudo-crossover function of the Sug- 
den factor which accounts for the MR description of 
the classical-to-critical crossover, in the homogeneous do- 
main (see the discussion in next section). The corre- 
sponding curves labelled Mi?+ in Figs. [1]: and[TJl, con- 
firm the perfect agreement with the master behavior of 
the one-component fiuid subclass when the asymptotic 
term of S (T*) corresponds to the one of 5*. (T*) for 
r* ^ 0. 
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(a) 


exponent 




S2 


i 






(b) 


exponent 






i 




Y 




V = 0.6303875 


2.121008 


22.9007 


1 


40.0606 


-0.098968 




7 = 1.2395935 


3.709601 


22.9007 


1 


29.1778 


-0.178403 




A = 0.50189 






2 


11.9321 


-0.15391 




A = 0.50189 






2 


11.7625 


-0.282241 




l^MF = 0.5 






3 


1.90235 


-0.00789505 




7AfF = 1.0 






3 


2.05948 


-0.0185424 














5.81623 












1},^ 


8.56347 



Table III: Values of the universal exponents and universal parameters for the crossover functions estimated in Ref. |25|]; Part 
(a) correlation length case in the homogeneous domain; Part (b) susceptibility case in the homogeneous domain. 



Correlation length 


Susceptibility 




f = 0.6303875 


7 = 1.2396935 






= 0.471474 


(Zj)~^ = 0.269571 


©{I/} ^ 4 288 10"^ 


Z+ = [Z+L<i« ( 


;ew>)y' = o.57 


z+ = [z+ (l^i^)"' (e<i^>)T' = 0.119 


^{1/} ^ 1.7410-4 


Z*^^'* = L"t^^> = 25.6988 


Z<^^> = (*(^^>)'1 " = 1950.7 


L*^-'^ = 25.6988 


A = 


0.50189 






Z^'+ = 0.68Z^'+ = 5.81623 


Zi'+ = 8.56347 






= 0.682:^'+ = 0.377 


^x+=^x^ (e<^-f>)''= 0.555 











Table IV: Universal and master constants of Egs . II32II a nd l|33p for the correlation length and the susceptibility, respectively, 
in the homogeneous domain (see text and Refs. [251. l35l| for details). Upper part (lines 1 to 4) refers to the Ising-like leading 
term; The values of the three characteristic numbers of the one component fluid "subclass" are reported in column 3, that 
demonstrates the unequivocal relation between the "master" crossover functions |3^ and the "MR" crossover functions (25l |. 
Lower part (lines 5 to 7) refers to the flrst term of the confluent correction to scaling. 



3.3. The master leading power law of the 
renormalized Sugden factor 

In the preasymptotic domain defined by Eq. ([39l) . 
the above formulation of the master singular behavior 
of S (T*), with T* > 0, can be approximated by a re- 
stricted (two term) expansion of equation 



1 



(41) 



where the decorated hat labels pseudo-physical quanti- 
ties. Equation l|4ip contains the asymptotic constraint of 
Eq. (|28|) . written following the master description 



lim 



Sir*) 



r'-.o+ 



= lim[5;.(|T*|)] 



(42) 



where 5*. (|r*|), with T* < 0, is given by Eq. while 
the difference occuring to the first order of the confiuent 
corrections to scaling is discussed below (see § 4.1). The 
leading amplitude Zg has the master form 



Zs = R 



(43) 



Using the universal values = 0.376 (±0.017) [7, 

[13, M, Hi, = 0.0574 (±0.0020) Rt = 



0.2696 (±0.0007) [23|, estimated for the Ising-like univer- 
sality class, and the values Z+ — 0.119, Z^ — 0.57 (see 
Table lEI, we obtain 



Zs = 2.47 (±0.17) 



(44) 



We note that the error-bar reported for each univer- 
sal amplitude combination only account for theoretical 
uncertainties on the estimated values of the universal 

while the "best" cen- 



combinations R^^ 



and i?^, 



tral values of the master amplitudes Z+ and Z^ are 
estimated using xenon as a standard one component 
fiuid . The large error-bar (±5%) on R^^ accounts 

0.367 (±0.009) and 



for the theoretical values R'^^ 



~ 0.372 (±0.009) estimated by Zinn and Fisher \M 



from numerical studies of three-dimensional Ising mod- 
els, the (min and max central) values i?+ ~ 0.36 (±0.01) 
and ~ 0.39 (±0.03) quoted by Privman et al 
on the basis of previous theoretical calculations, and 
the median values ~ 0.386 (±0.1) and R+^ ~ 
0.381 (±0.01) [l3| which were initialy obtained from the 
analysis of the experimental situation for fluids (see Refs. 
[ll],a3,Iii, 15, 16, 17, 18, 19, 20, 21, 2I). 

The published data of the effective exponent-amplitude 
pair {(pe',So,e} reported in Table U (colums 3 & 4) allows 
one to validate this leading master description at finite 
distance to the critical point, using a method equiva- 
lent to the one proposed by Moldover [l3| to estimate 
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So by averaging the values of ^^,|0.3-is in the vicinity 
of |Ar*| = 0.01. The corresponding Moldover's val- 
ues (noted Sq^^ to recall for the use of the theoretical 
value Lp = 0.935), are given in column 5 of Table [H In 
our present work, we have estimated So^ip by the follow- 
ing relation So,ip = So^e (0.01)'''''" ' (see also column 
5, Table [H. From these "measured" amplitude data at 
I At* I = 0.01, the corresponding calculated values (col- 
umn 6) of Zs.^ = {ac)^~'^ (g*)^ (^c)~^ (5^c)~'^ -5*0,^ [see 
Eq. (|22|) ]. are in close agreement with the asymptotic 
limit = 2.47 estimated from above hyperscahng con- 
siderations. The mean value of the data reported in col- 
umn 6 is (Zs^tp) = 2.450. The residuals SZg^tp (column 
7), expressed in %, are of the same order of magnitude 
(±3.1'%) than the experimental uncertainty (±5'%) [see 
for example the review of Moldover [l3[ for a detailed 
analysis of the realistic experimental errors]. 

This extended master behavior is illustrated in Figs, 
[ifc and [Hi by the curve labelled MRq which corresponds 
to the pure power law of equation 



4.1. 



4. DISCUSSION 

Ising-like universal features within the 
preasymptotic domain 



As demonstrated in Refs. |25l.l35|. each crossover func- 
tion obtained from the MR scheme can be approximated 
by a restricted (two term) Wegner-like expansion in the 
Ising-like preasymptotic domain which extends up to 



IT* I < C 



{1/} 

PAD 



/rising 
'-'PAD 



510" 



(see the corresponding arrows in |T*| axis of Fig. [TJ; 
and [Hi). Therefore, in addition to Eq. l(2T|) related to 
the master singular behavior of the renormalized Sug- 
den factor, we are also interested by the following similar 
equations 



M 



LV 



-■M 



IT* 



1 



(46) 



Zs \T 



* \v 



(45) 



where Zs = 2.47 [see Eq. (|44])]. In Fig. [TJi, the two 
lines labelled up [Eq. l(44|) with Zs = 2.64] and down 
[Eq. (|44l) with Zs — 2.30], respectively, account for the 
theoretical error-bar attached to this central value of Zs- 
Therefore, at least for a temperature-like range such that 
|T*| < 0.1, all the experimental results measured at fi- 
nite temperature distance to the critical point appears 
"condensed" within these two lines. As noted previously, 
such a good agreement result from the "universal" me- 
dian value ife = ifi — 0.935 of the effective exponent in 
the vicinity of At* — 0.01. De facto, the asymptotical 
universal features can be observed in an extended asymp- 
totic domain, since the confluent corrections to scaHng 
attached to the exponent A are i) only governed by the 
single scale factor Yc whatever the singular property (as 
already shown for the correlation length, the susceptibil- 
ity, and the order parameter density), and, ii) certainly 
very small in amplitude for the Sudgen factor case. How- 
ever, the present theoretical and experimental levels of 
uncertainties are of same order of magnitude and remain 
too high to provide an accurate estimation of the sign 
and the amplitude of these (small) confluent corrections. 

As our explicit Eq. I|40p is restricted only to the uni- 
versal features related to hyperscahng, there is a need for 
theoretical studies in the future to directly estimate the 
classical-to-critical crossover of the surface tension and 
the Sudgen factor in the non-homogeneous domain. An- 
ticipating these investigations, the following discussion 
gives some complementary quantitative evaluations on 
the extended temperature-like range where the asymp- 
totic leading power law of Eq. l(45|) can be correctly used 
to predict the Sudgen factor behavior (since the applica- 
bility of the scale dilatation method goes far beyond the 
one of the unvalid corresponding state principle) . 



E* = Zy^ IT* 



Zl\T* 



(47) 



related to the master singular behaviors of the renor- 
mahzed order parameter density [see Eq. lfT6| ] and 
renormalized surface tension [see Eq. ifTT)) ]. respectively. 
Obviously, the hyperscahng law dv — j + 2f3 provides 



the universal combination 



(Zm) 



— Rr^ ( R 



while Eq. i|19p provides the "trivial" relation Zs = 
Both of these amplitude combinations relate unequiv- 
ocally Zm and Zy to the selected characteristic lead- 
ing amplitudes Z+ and Z^ of the one-component fluid 

subclass. Alternatively, Zs and Z^ are unequivocally 

related by the universal amplitude combination = 

Zs (2^) ■ In such a case, we can also calculate 



the universal values Zs = R 



d-1 



= 1.750 and 



Zs = 2^ = 1.867 of the corresponding leading ampli- 
tudes for the respective crossover functions estimated in 
the MR scheme [with Z+ = 2.121, Z+ = 3.7096, and 



Za/ — 



= 0.9375; see Ref. M for 



detail]. Furthermore, in the relations [similar to Eqs. 
((32l) and (|33l) ] which deflne the master crossover func- 
tions for the order parameter density, the surface tension 

1 1 f I 

and the Sugden factor, the respective prefactors Z]^/^, 
Z^ , and Z^ account for their unequivocal estima- 
tion only using the three characteristic numbers L^^-''^, 
e^i^>, and ^'{i^^ of the one-component fluid subclass, 
such that, 



Z 



{1/} _ 

M ^ 



Zj\4 



(48) 
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Order parameter density 


Interfacial tension 


Sugden factor 


P — U.oZ0/o40 


J. o., 1 Ofin77c: 

— zt^ — l.ZDU/ (0 


= 


^<p- 13 = 0.9349905 


= (-RJZJ) ( 1 = 0.937528 


Ze = -RJ"^ (Z^ j = 1.6915 


Zs = i?± 




' = 1.8042 




Zm = Zm (l^i^)" v1>{i« fe<^«V = 0.468 
V / V / 


2e=Ze fL{i«V"7eW}V = 1.1558 


Zs^ 




T^il/} /'tT .T/{1/} 9Q/IS7S 


= ( IL"^ M = DDU.4ZO 






= 223.634 


(z(">)° 








— 1 










A = 0.50189 






Zli = 0.9Z^'+ = 7.70712 






Zlt = Zjv, (e<i^>)'^ = 0.9^^'+ ^ 0.5 


2^^^!, ^#-1 (seeSg.l) 
1 

4=F ~ 0.9 -^Z^^ 0.5 





^0 

0.9^ ^0 


(see -Fi*;. 1) 



Table V: Universal and master constants for the order parameter density (column 1), the surface tension (column 2), and the 
Sugden factor (column 3), in the nonhomogeneous domain. Upper part (lines 2 to 6) refers to the Ising-like leading term; The 
unity value of the combinations between the master prefactors reported in line 5 demonstrates that the asymptotic master 
crossover agrees with the two-scale-factor universality of the Ising-like systems. Lower part (lines 7 to 9) refers to the first term 
of the confluent correction to scaling (see text for detail). 



z. 



{1/} 



Z+ (9(1/}) 



Z 



{1/} 



Zs 



zs(e{i 



(49) 



(50) 



Equations l(48|) to (|50l) close the master representation 
of the singular behavior of the renormalized interfacial 
properties in the nonhomogeneous domain, in agreement 
with the two-scale factor universality of the Ising-like sys- 
tems [see the corresponding values of the universal and 
master quantities reported in Table |Vj. 

Now, using Eq. (flOl) to compare the respective first 
confluent amplitudes of Eqs. (061), we obtain 

Zg = — Z\j. >From Fig. [TJl, the asymptotic master 
singular behavior expected for S*, {\T*\) is compatible 
with the following universal values of the corresponding 
amplitude ratios 



2l7 



1 



4f =0.9, 







(51) 



Such hypothesized "universal ratios" of Eq. (jSTl) are con- 
sistent with Ising-like universal features of the asymp- 
totic crossover estimated from the MR scheme, which are 
only characterized by a single confluent amplitude within 
the Ising-Hke preasymptotic domain. Here these univer- 
sal features are preserved via the universal ratio value 

~ 0.9, selecting -2?;'+ as a characteristic confluent 

amplitude (see § 3.2 above). However, it is also impor- 
tant to note that this expected crossover must satisfy 



the scaling law (/3 = — in the infinite limit |T* | oo, 
which leads to the mean field value (fiMF = 1, using the 
mean-field values /3mf = 5 and c/jaif = § [4a|. In the 
range |T*| > 00, the experimental results reported in 
Fig. [1^ to [Hi are in disagreement with such a mean-field 
prediction (see also below the § 4.3). 

In addition, we note that the hyperscaling descrip- 
tion using a pseudo-crossver function issued from sin- 
gular properties in the homogeneous domain generates 
uncorrect results in the complete temperature range, i.e. 
from the first-order contribution of Ising-like confiuent 
exponent A until the leading contribution related to the 
mean-field exponent lpmf- 

For example, in our scheme based on the hyperscal- 
ing law tp = [see Eq. l(29|) ]. the confiuent ampli- 

in Eq. (|4T1) can be made equal to 

'■ 0.466, leading to a universal ratio 



tude Z 



s 

Zb+ - 



4" 



= i 1 



0.84 which is different from zero. 



Similarly, a description based only on the hyperscaling 
law if = 2v — (i [see Eq. l(24l) ] needs to replace the inter- 
facial tension by the inverse squared correlation length 
in Eq. (fT9|) , and provides another pseudo-crossover func- 
tion, given by the equation 



5(|T1) = <x- 



1 



LV 



(|r*|) 



1 



n 2 



t (T*) 



(52) 



where the decorated tilde distinguishs new pseudo- 
physical quantities from those of Eq. (f40|) . In that case, 
a mixing occurs between properties in the homogeneous 
(£* (T*)) and nonhomogeneous (A^^y(|T*|)) domains. 
In the Ising-like preasymptotic domain, accounting for 



the relation with Zs — R^^Zm {^Z'^^ , Eq. 



(1521) can 
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be approximated by 

S{T*) = Zs{\T*\f 



Zl{\T* 



(53) 



In this latter scheme, the confluent ampUtude Zg in Eq. 
([53l) was estimated equal to Zl = Z\,r + 2Z^+ w 1.254, 



^1 

leading to a universal ratio --^ 



= 0.9 



2.26 



which is also signiflcantly different from zero. 

Looking now to the contribution of the leading term 
close to the Gaussian fixed point, our pseudo-crossover 
functions estimated above does not account for the ap- 
propriate mean-field-like description due to the failure 
of the two hyperscahng laws ip = (which gives un- 
correct value <fMF — f ) and ip = 2v — (3 (which gives 
uncorrect value ipuF ~ 5) when we use the correspond- 
ing mean-field values 7mf = 1, vmf — \ and fihiF — \ 



extended asymptotic domain (labelled EAD) goes up to 
the limit 



IT* I < « 0.03 



(55) 



(see the corresponding arrow noted C^^^^ in |T*| axis). 
Within |T*| < L^^^^ ^ the observed master behavior can 
be well-represented by 5*. = Zs \T*f [see Eq. (05])], 
in conformity with the Ising-like universal features esti- 
mated from the MR scheme. We note that such Ising-Hke 
nature of S*, in this extended |T*| range complements 
in a self-consistent manner our previous analysis [W] of 
the master behavior of the renormalized order parameter 
density along the VLE fine. 



4.2. Ising-like master behavior in the extended 
asymptotic domain 

In spite of the absence of accurate theoretical mod- 
elling for interfacial tension and Sugden factor along the 
VLE fine, the MR description of the master crossover ob- 
served for the one-component fluid subclass can be used 
to provides a reasonable estimation of the renormalized 
correlation length in the nonhomogeneous domain, using 
the following equation 

t (r*<0) = ^* (54) 

where t (T* > 0) of Eq. ([32]) is the renormalized cor- 
relation length in the homogeneous domain. Eq. l(54|) 

assumes that the universal ratio ^^l^.^gj = 1-96 is inde- 
pendent of the renormalized temperature like field. The 
result (for T* < 0) is illustrated as a £* (T* < 0) gradua- 
tion of the upper horizontal axis of Figs, [ifc and [Hi. We 
recall that I* gives the best estimate of the ratio be- 

O etc 

tween the effective size (^) of the critical fluctuations and 
the effective size (etc) of the attractive molecular interac- 
tion, the latter one being approximated by the dispersion 
forces in Lennard-Jones-Hke fluids which extend over a 
short range slighly greater than twice the equilibrium dis- 
tance Te between two interacting particles of flnite hard 
core size cr (thus etc ~ 2re, with re ^ cr). Therefore, 
(|T*| = Ccic) ~ 1 in the upper axis of Figs[lfc and [Hi 
is a rough estimate of the microscopic range of the molec- 
ular attractive interaction between fluid particles. Such a 
thermal field limit corresponds to the value Ccic ~ 0.15 
(here the supscript CIC recall that the effective extend 
of the short-ranged molecular interaction corresponds to 
the size of the critical interaction cell) . Looking then to 
the "Ising-Hke" nature of S*, (|T*|), we observe in Figs, 
[ifc and [TH a noticeable extension of the critical range as- 
sociated to the condition (.* {\T*\) > 3. Therefore, the 



4.3. Non-critical behavior beyond the Ising-like 
extended asymptotic domain 

In Ref. [H3|, it was observed for the xenon case, 
that the real crossover for the effective exponent Pe ap- 
pears in the thermal field range \T>qq\ w 0.1 — 1 where 
£* {\T*\) = 1 (see Fig. [TJl). In terms of comparison be- 
tween the correlation length and the range of the micro- 
scopic intermolecular interaction, the situation is similar 
to the one encountered in the homogeneous domain for 
the real crossover for the effective exponent 7e [s^l- In- 
deed, when £* < 1, any MR crossover function is not 
appropriate to account for the real non-universal behav- 
ior of the one-component fiuids. We recall for example 
that £* = ^ (see the limiting curve m in Fig. [TJi) corre- 
sponds to a (non-master) microscopic arrangement where 
the direct correlation distance between two interacting 
particles is « Ve (i.e., £ (At* < 0) « re ^ ct). As pre- 
viously noted in Ref. [57||, the nonhomogeneous fluid is 
then made of coexisting gas and liquid which show sig- 
nificant differences in the averaged quantity of particles 
inside the critical interaction cell. Moreover, these differ- 
ences increase approaching the triple point temperature, 
since the low density gas tends to behave as a perfect gas 
with one (i.e. non-interacting) particle within the CIC 
volume, while the condensed liquid tends to minimize the 
configuration energy of one particle by enclosing them in 
a particle cage made with an increasing number (up to 
twelve for rare gas case) of the closest neighboring (re- 
pulsive) particles (i.e. the mean size d of the cage is such 
that Tg < d ~ cr). For such "low" and "high" local den- 
sities, cooperative density fiuctuations have no physical 
sense at length scale larger than ac and the non-universal 
characteristics of each fluid are only involved in the ther- 
modynamic properties, as clearly illustrated in Fig. [T}! 
for the Sugden factor case by the significative increasing 
differences between the rescaled data for xenon and water 
in the range \T*\ > 0.2. 
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5. CONCLUSION 

We have provided an asymptotic description of the sin- 
gular behavior of the renormaUzed Sugden factor (i.e. 
the renormalized squared capillary length) of the one- 
component fluid subclass. This master crossover behav- 



ior can be observed up to |T*| « 0.03 (or i* ~ 3) in 
the non-homogeneous domain, as already noted for the 
renormalized order parameter density. In a future work, 
we will show that this master crossover behavior can be 
useful to estimate the parachor correlation along the VLE 
line. 
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